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In  inertial  confinement  fusion  (ICF),  the  possibility  of  ignition  or  high  energy 
gain  is  largely  determined  by  our  ability  to  control  the  Rayleigh-Taylor  (RT)  instability 
growth  in  the  target.  The  exponentially  amplified  RT  perturbation  eigenmodes  are  formed 
from  all  sources  of  the  target  and  radiation  non-uniformity  in  a  process  called  seeding. 
This  process  involves  a  variety  of  physical  mechanisms  which  are  somewhat  similar  to 
the  classical  Richtmyer-Meshkov  (RM)  instability  (in  particular,  most  of  them  are  active 
in  the  absence  of  acceleration),  but  differ  from  it  in  many  ways.  In  the  last  decade,  the 
radiographic  diagnostic  techniques  have  been  developed  that  made  direct  observations  of 
the  RM-type  effects  in  the  ICF-relevant  conditions  possible.  New  experiments  stimulated 
the  advancement  of  the  theory  of  the  RM-type  processes.  The  progress  in  the 
experimental  and  theoretical  studies  of  such  phenomena  as  classical  RM  instability  in 
finite-thickness  targets,  re-shock  and  re-rarefaction  of  the  RM-unstable  material 
interfaces,  ablative  RM  instability,  feedout,  and  perturbation  development  associated  with 
impulsive  loading  is  reviewed. 

Keywords:  inertial  confinement  fusion,  hydrodynamic  instability,  Richtmyer-Meshkov 
instability,  shock  wave,  expansion  wave 

1.  Introduction 

Interfacial  instabilities,  such  as  Rayleigh-Taylor  (RT)  (Rayleigh  1883,  Taylor 
1950)  and  Richtmyer-Meshkov  (RM)  (Richtmyer  1960,  Meshkov  1969),  and  turbulent 
mixing  driven  by  these  instabilities  are  ubiquitous  in  nature.  Their  manifestations  range 
from  a  boiling  tea  kettle  in  the  kitchen  to  Gulf  Stream  in  the  ocean  to  exploding 
supernovae  (Schmidt  2006).  Studies  of  instabilities  and  mixing  have  become  one  of  the 
most  actively  developed  areas  of  fluid  dynamics.  Rapid  advancement  in  flow  diagnostics 
provided  a  vast  amount  of  data  from  numerous  experiments  with  liquids  and  gases  on 
channel,  falling-  or  accelerated-tank  facilities  and  shock  tubes,  see  Niederhaus  &  Jacobs 
2003,  Holder  et  al.  2003,  Sadot  et  al.  2005,  Chapman  &  Jacobs  2006,  Mueschke  et  al. 
2006,  Dimonte  et  al.  2007,  Motl  et  al.  2007,  Orlov  &  Abarzhi  2007,  Kumar  et  al.  2007 , 
Schwaederle  et  al.  2007,  Wilkinson  &  Jacobs  2007,  Scase  et  al.  2008,  and  references 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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therein.  Progress  in  massively  parallel  numerical  simulations  (Cabot  &  Cook  2006),  as 
well  as  advancements  in  the  theory  greatly  improved  the  general  understanding  of 
mixing.  The  physical  mechanisms  responsible  for  the  perturbation  growth  and  turbulent 
mixing  in  shock  tubes,  channel  facilities  and  imploding  laser  fusion  pellets  are  similar  in 
many  ways,  which  is  why  the  knowledge  gained  in  the  room-temperature  experiments 
proved  helpful  for  the  inertial  confinement  fusion  (ICF),  guiding  the  target  design  and 
facilitating  the  data  analysis. 

At  the  same  time,  understanding  of  many  issues  related  to  the  perturbation  growth 
and  mixing  in  the  ICF  targets  requires  experimental,  theoretical  and  computational 
studies  performed  specifically  for  the  ICF-relevant  conditions,  in  high-energy-density 
(HED)  plasmas  at  or  above  solid-state  density,  heated  to  temperatures  from  eV  to  keV 
and  compressed  to  multi-Mbar  pressures.  Implosions  of  ICF  capsules  include  processes 
and  regimes,  which  are  hard  or  even  impossible  to  emulate  in  a  different  environment, 
such  as  complex  history  of  acceleration  and  deceleration  of  its  shell(s)  and  interface(s), 
their  multiple  re-shocks  and  re-rarefactions  (Lindl  1998,  Bodner  et  al.  1998,  Atzeni  & 
Meyer- ter- Vehn  2004,  Pfalzner  2006).  Much  of  the  perturbation  seeding  and  growth  in 
such  targets  occur  at  the  ablation  fronts,  which  can  only  be  produced  by  interaction  of 
intense  radiation  or  hot  plasma  with  solid  matter  and  whose  properties  are  very  different 
from  those  of  material  interfaces. 

Although  the  HED  plasmas  are  a  challenging  environment  for  observing 
hydrodynamic  instabilities,  some  of  the  relevant  basic  hydro  phenomena,  particularly  the 
RM  instability  growth  and  related  effects,  are  actually  easier  to  observe  and  analyze  in 
the  ICF/HED  plasma  experiments.  Here  are  some  examples. 

1)  There  is  a  clear  advantage  in  experimenting  with  clean,  well-characterized 
rippled  interfaces  separating  two  different  solid  materials.  A  solid  target  can  remain  at 
rest  until  shocked  (Dimonte  &  Remington  1993,  Dimonte  et  al.  1996,  Glendinning  et  al. 
2003,  Aglitskiy  et  al.  2006).  In  conventional  shock-tube  experiments  with  two  gases,  the 
presence  of  thin  membrane  affects  the  early-time  perturbation  growth.  Since  the  original 
work  of  Meshkov  1969  up  to  the  most  recent  publications  (Erez  et  al.  2000,  Mariani  et  al. 
2008)  this  “membrane  effect”  proved  to  be  difficult  to  eliminate  or  accurately  account 
for. 

2)  The  rarefaction-triggered  counterpart  of  the  RM  instability,  developing  when  a 
rarefaction  wave  reaches  a  material  interface  from  a  high-density  side,  has  been 
theoretically  discovered  long  ago  (Li  &  Book  1991,  Li  et  al.  1991)  and  is  undoubtedly 
important  for  ICF.  No  room-temperature  facility  is  available  for  its  experimental  study.  In 
an  ICF  experiment,  launching  a  rarefaction  wave  is  easy,  it  happens  naturally  when  a 
radiation-driven  shock  wave  breaks  out  at  the  rear  surface  of  the  target,  as  in  Smitherman 
et  al.  1999,  Shigemori  et  al.  2000,  Aglitskiy  et  al.  2001b,  2006. 

3)  Most  important  manifestation  of  the  perturbation  growth  in  a  radially  imploded 
ICF  target  is  the  lateral  mass  redistribution  inside  it  that  reduces  the  target  uniformity  and 
thus  degrades  the  fusion  neutron  yield.  This  process,  which  is  hard  to  study  in  room- 
temperature  experiments,  is  routinely  observed  in  the  ICF-relevant  conditions,  which  is 
made  possible  by  advanced  x-ray  radiographic  diagnostic  techniques. 

The  RM-type  hydrodynamic  processes,  besides  the  classical  RM  instability, 
include  a  variety  of  related  phenomena,  some  of  which  are  specific  to  the  ICF  targets, 
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while  the  other  are  of  more  general  nature.  We  focus  on  the  basic  physical  mechanisms 
responsible  for  the  onset  and  early-time,  small-amplitude  perturbation  evolution.  This  is 
because  the  interest  to  the  RM-type  processes  in  the  ICF  is  determined  by  their  role  in  the 
RT  “seeding,”  i.  e.,  formation  of  the  RT  eigenmodes,  which  are  subsequently  amplified 
by  the  exponential  RT  growth,  from  the  initial  target  imperfections  and  radiation  non¬ 
uniformities.  No  effort  is  spared  to  make  the  ICF  targets  and  the  radiation  fields  as 
smooth  as  possible,  so  the  initial  amplitudes  of  all  the  Fourier  modes  representing  their 
non-uniformity  are  very  small.  Many  of  the  RM-type  processes  are  oscillatory  and  can 
cause  only  limited  growth  of  initial  perturbation  growth.  The  other,  like  the  classical  RM 
instability,  can  cause  large  growth,  but  the  ICF  target  designs  and  implosion  scenarios 
typically  ensure  that  any  RM-type  growth  is  limited  to  one  order  of  magnitude  or  so 
before  the  RT  growth  begins. 

For  this  reason,  the  small-amplitude  RT  seeding  is  hard  to  observe  directly  and 
notoriously  difficult  to  model  numerically,  see  Zalesak  et  al.  2005.  Until  very  recently, 
little  experimental  data  were  available  for  testing  the  codes  and  theories  in  the  RT- 
seeding  relevant  regimes.  The  situation  changed  with  the  development  of  new  diagnostic 
techniques  (Aglitskiy  et  al.  2001a,  2001b,  2002,  2006,  Gotchev  et  al.  2006),  which  made 
it  possible  to  study  the  whole  array  of  the  previously  unexplored  RM-type  processes  in 
the  ICF  targets. 

In  this  article  we  summarize  the  results  of  such  studies  performed  over  the  last 
decade,  mainly  at  the  Naval  Research  Laboratory  (NRL)  and  based  on  direct 
observations,  most  of  which  have  been  done  for  the  first  time.  In  Section  2  we  review  the 
drivers  and  diagnostic  capabilities  available  for  ICF-HEDP  experiments  on  RM-type 
instabilities.  In  Section  3  we  discuss  the  so-called  ablative  RM  instability  which 
approximates  the  limiting  case  of  Atwood  number  A  =  1  of  the  classical  light-to-heavy 
RM  instability,  when  a  rippled  ablation  front  is  simultaneously  produced  and  shocked  by 
a  laser  or  x-ray  radiation  rapidly  turned  on.  Section  4  is  dedicated  to  the  opposite  limiting 
case  of  the  classical  heavy-to-light  RM  instability  A  =  -1 ,  the  feedout  situation,  when  the 
shock  wave  arrives  to  a  rippled  rear  surface  of  a  target,  behind  which  there  is  no  low- 
density  material.  In  Section  5  we  describe  the  effect  of  a  re-shock,  which  stops  the 
acceleration  of  a  RT  unstable  ablation  front  and  quenches  the  RT  growth.  In  Section  6  we 
discuss  the  perturbation  evolution  in  a  shock-rarefaction  flow  produced  by  impulsive 
loading:  A  short  powerful  laser  pulse  that  is  turned  off  immediately  after  depositing  a 
finite  energy  on  a  rippled  target  surface.  In  Section  7  we  conclude  with  a  discussion. 

2.  Laser  drivers  and  diagnostics 

The  ICF-relevant  high  energy  density  parameter  range  for  hydrodynamic 
experiments  means  multi-Mbar  pressures  (Lindl  1998,  Atzeni  &  Meyer- ter- Vehn  2004, 
Pfalzner  2006.)  This  is  the  pressure  range  required  for  strong  shock  compression  of  solid 
targets,  see  Zel’dovich  &  Raizer  2002,  Chapter  XI.  The  multi-Mbar  pressure  needs  to  be 
maintained  for  a  sufficiently  long  time,  over  a  sufficiently  wide  area  to  produce  the 
required  ID  flow  and  observe  the  perturbation  development  characteristic  of  it.  These 
conditions  determine  the  requirements  on  laser  drivers  suitable  for  the  relevant 
experiments. 
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The  ablative  pressure  produced  by  direct  laser  irradiation  of  a  solid  deuterium  or 
plastic  target  can  be  estimated  as 


f  T  \ 
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AlJ 


2/V  A 
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V  re  J 
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“1 2/3 


Mbar, 
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where  I  and  AL  are  the  laser  intensity  and  wavelength,  respectively,  mu  is  the  atomic 
mass  unit  and  re  is  the  classical  electron  radius,  see  Manheimer  et  al.  1982.  For  indirect 
soft  x-ray  drive,  the  ablative  pressure  in  the  hohlraum  cavity  is  related  to  its  radiation 
temperature  Tr  by 

p  m  3  x  1CT7  x  [Tr  (eV)]7/2  Mbar,  (2) 

see  Lindl  1998.  Laser-produced  pressure  of  p  =  10  Mbar  corresponds  to  the  intensity 
/ranging  from  31  TW/cm  to  66  TW/cm  for  laser  wavelength  varied  from  0.248  pm 
(KrF  laser)  to  0.53  pm  (second  harmonic  of  Nd-glass  laser  radiation).  The  hohlraum 
temperature  required  to  produce  the  same  pressure  is  140  eV. 

Table  I  presents  the  parameter  range  for  the  RM-type  hydrodynamic  experiments 
performed  so  far  in  the  ICF-relevant  conditions.  It  is  relatively  narrow:  driving  pressure  p 
from  7  to  45  Mbar,  laser  pulse  duration  tL  from  1.5  to  11  ns,  shock  velocity  D  from  20  to 
60  pm/ns,  payload  diameter  from  300  to  800  pm,  ripple  wavelength  A  from  20  to  150 
pm.  Indeed,  estimating D  s  ^ p/ p  ,  we  find  that  a  ~10  Mbar  pressure  drives  a  shock 
wave  into  a  solid  plastic  or  a  light  metal  ablator  like  Be  (Dimonte  &  Remington  1993)  or 
Al  (Smitherman  et  al.  1999)  at  the  velocity  ~  2-3xl06  =20-30  pm/ns.  To  make  the 
RM-type  phenomena  observable,  the  driving  pressure  in  most  cases  needs  to  be 
maintained  while  the  shock  propagates  more  than  one  ripple  wavelength  A .  This  imposes 
a  lower  limit  on  the  duration  of  the  pressure  drive.  On  the  other  hand,  the  accelerated  area 
should  be  wide  enough  to  cover  a  large  number  of  wavelengths  A ,  to  ensure  periodicity. 
These  two  conditions  require  high  laser  energy  delivered  to  the  target  on  a  ns  time  scale. 
In  direct-drive  experiments,  like  those  reported  by  Aglitskiy  et  al.,  2001a,  2001b,  2002, 
2006,  Glendinning  et  al.  2003,  Weaver  et  al.,  2004,  Gotchev  et  al.,  2006,  about  0.5-2. 5  kJ 
of  laser  energy  needs  to  be  deposited  into  the  focal  spot.  For  soft  x-ray  drive  experiments, 
laser  energy  of  20-30  kJ  has  to  be  delivered  into  the  hohlraum,  because  the  ablator 
pushing  the  payload  can  constitute  only  a  small  fraction  of  its  wall  area  (Dimonte  & 
Reminton  1993,  Dimonte  et  al.  1996). 

Only  high-power,  high  energy,  multi-beam  laser  systems  developed  for  the  ICF 
can  satisfy  these  requirements.  The  10-beam,  40-kJ  NOVA  laser  at  the  Lawrence 
Livermore  National  Laboratory  (LLNL)  had  been  shut  down  in  1999,  but  all  the  other 
drivers  listed  in  Table  I  remain  in  operation:  the  12-beam  GEKKO-XII/HIPER  laser 
facility  at  the  Institute  of  Laser  Engineering,  Osaka  University,  Japan  (2.4  kJ  in  0.35 1  pm 
wavelength  radiation,  see  Yamanaka  1985),  the  60-beam  OMEGA  laser  facility  at  the 
Laboratory  for  Laser  Energetics,  University  of  Rochester  (30  kJ  in  0.351  pm,  see  Boehly 
et  al.  1995)  and  the  56-beam  Nike  KrF  laser  facility  at  the  Naval  Research  Laboratory  (3 
kJ  in  0.248  pm,  see  Obenschain  et  al.  1996). 
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Next-generation  high  energy  lasers  are  presently  being  built,  opening  new 
opportunities  for  hydrodynamic  experiments.  Construction  of  the  196-beam,  1.8  MJ  in 
0.351  pm  National  Ignition  Facility  (NIF)  at  LLNL  has  been  completed  in  March  2009. 
A  similar  ignition-scale  facility  Laser  Megajoule  (LMJ,  240  beams,  1.8  MJ  in  0.351  pm) 
is  under  construction  at  CESTA,  France  (Besnard  2007).  A  4-beam,  10  kJ  in  0.351  pm 
LMJ  prototype,  Ligne  d’lntegration  Laser  operates  in  France  since  2002  (Eyharts  et  al. 
2006).  A  8-beam,  15-20  kJ  in  0.351  pm  glass  laser  system  operates  in  China  since  2005 
as  a  prototype  of  the  64-beam,  150-200  kJ  Shenguang-III  facility  (He  &  Zhang  2007).  A 
4-beam,  9  kJ  in  0.351  pm  glass  laser  system  Luch  is  under  construction  at  VNIIEF, 
Russia  (one  of  its  beams  is  operational  since  2002,  see  Garanin  et  al.  2005)  as  a  prototype 
of  the  128-beam,  300  kJ  Iskra-6  facility. 

Table  I.  Parameters  of  the  RM-type  hydrodynamic  experiments  on  high  energy  laser  facilities. 
Here  AL  and  tL  are  laser  wavelength  and  pulse  duration,  respectively,  p  is  the  pressure 
maintained  by  the  driver,  0  is  the  diameter  of  the  laser-  or  soft-x-ray  driven  area,  D  is  the 
velocity  of  the  shock  wave  driven  into  the  target,  A  is  the  ripple  wavelength,  hv  is  the  photons 
energy  of  backlighting  x  rays,  with  backlighter  material  indicated. 
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The  RM-type  hydrodynamic  effects  involve  redistribution  of  mass  in  cold  and 
dense  areas  of  laser-driven  targets.  Such  mass  motions  are  observed  directly  using 
absorption  of  hard  x  rays  as  a  diagnostic  tool.  For  this,  a  powerful  external  source  of  x 
rays  is  needed  to  illuminate  the  target,  and  an  imaging  system  -  to  record  the  distribution 
of  mass  probed  by  the  x  rays.  The  ICF  x-ray  sources  are  backlighters  -  small  targets 
made  of  high-atomic-number  materials,  which  are  irradiated  by  the  separate  laser  driver’s 
beams  to  convert  laser  energy  into  keV  x-ray  radiation.  Table  I  lists  the  backlighter 
materials  and  the  corresponding  energies  of  the  x-ray  photons  used  in  the  RM-type 
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experiments.  To  ensure  that  the  hard  x-ray  source  is  bright  enough  to  probe  thick,  dense 
targets,  experiments  often  require  a  large  amount  of  laser  energy  for  backlighting  -  about 
as  much  as  is  used  to  drive  the  target  (Gotchev  et  al.  2006)  or  even  more  (Dimonte  & 
Remington  1993,  Dimonte  et  al.  1996,  Glendinning  et  al.  2003). 

For  any  backlighting  scheme  there  are  two  desirable  features  of  the  imaging 
system:  smaller  size  of  the  backlighter  source  and  filtration  of  the  imaging  radiation  from 
the  self-emission  of  the  object.  If  the  dimension  of  the  source  of  backlighting  radiation  is 
small,  the  laser  power  required  to  create  it  can  be  much  smaller  than  the  power  driving 
the  object.  If  spectral  filtration  is  available,  the  surface  brightness  of  backlighter  must  be 
larger  than  the  object  brightness  only  in  the  bandwidth  of  the  spectral  line  being  used  to 
image  the  object.  The  large  spot  size  of  the  main  driver  is  necessary  to  minimize  the  edge 
effects  from  the  acceleration  process  and  therefore  requires  a  diagnostic  field  of  view  of 
approximately  1  mm,  so  that  the  entire  focal  spot  as  well  as  the  remaining  intact  part  of 
the  target  foil  are  within  the  picture. 

A  monochromatic  x-ray  imaging  system  based  on  Bragg  reflection  from 
spherically  curved  crystals  (Brown  et  al.  1997a)  has  been  successfully  used  on  the  NRL 
KrF  Nike  laser  to  observe  the  evolution  of  target  mass  perturbations.  Compared  to  other 
x-ray  imaging  techniques  like  pinholes,  Kirkpatrick-Baez  microscopes,  and  Fresnel 
lenses,  a  curved  crystal  imager  has  several  advantages.  It  is  a  high  throughput,  high 
spatial  resolution,  and,  by  nature,  monochromatic  diagnostic  technique  with  a  large  field 
of  view,  all  at  the  same  time. 
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Figure  1.  The  general  scheme  of  the  monochromatic  x-ray  imaging  diagnostics  fielded  on  Nike 
laser.  Also  shown:  the  shape  of  a  Si  backlighter  spectral  line  and  test  images  obtained  on  x-ray 
film  (magnified  in  the  inset). 

The  optical  scheme  of  monochromatic  x-ray  backlighting  (Rode  et  al.  1990, 
Pikuz,  S.  A.  et  al.  1995a,  1995b,  1997,  Pikuz,  T.  A.  et  al.  1995,  Aglitskiy  et  al.  1996, 
Brown  et  al.  1997a,  1997b,  Sanchez  del  Rio  et  al.  1997,  Pawley  et  al.  1997)  is  presented 
in  Fig.  1.  Imaging  of  just  top  and  bottom  points  of  the  main  target  is  shown.  For  the 
scheme  description  we  use  here  the  well  known  Rowland  circle,  which  is  a  curve  where 
all  x  rays  with  a  given  wavelength  A  focus  to  a  point.  There  is  no  restriction  on  the 
source  location  provided  the  radiation  wavelength  A  satisfies  the  Bragg  equation 
2d  sin#  =  in  A ,  where  d  is  the  crystal  interplanar  distance,  m  -  the  order  of  reflection,  and 
0-  the  Bragg  (grazing)  angle. 
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The  x-ray  backlighter  source  is  located  between  the  main  target  and  the  Rowland 
circle  of  radius  RI2,  where  R is  the  curvature  radius  of  the  crystal.  When  a  small  x-ray 
source  backlights  a  specific  point  of  the  object,  the  area  of  the  crystal  involved  in  the 
imaging  of  this  point  is  much  smaller  (0.5-1  mm)  than  the  crystal  itself.  Due  to  the 
focusing  of  the  backlighter  radiation  source  between  the  crystal  and  the  object  image,  it  is 
possible  to  protect  the  detector  from  a  significant  portion  of  the  self-radiation  of  the 
object  being  imaged.  Initially,  the  curved  crystal  imager  has  been  used  at  NRL  in 
combination  with  a  framing  camera  to  provide  sequences  of  high  spatial  resolution  still 
images  (up  to  four  snapshots  that  correspond  to  four  different  times)  taken  with  a  time 
resolution  of  200  ps  (Brown  et  al.  1997b).  This  diagnostic  technique  is  fully  adequate 
when  the  observed  mass  non-uniformity  varies  monotonically  in  time,  e.  g.,  during  the 
fast  RT  growth  (Brown  et  al.  1997b,  Sanchez  del  Rio  et  al.  1997). 


Observation  of  essentially  non-monotonic  evolution  of  the  processes  under  study 
like  the  oscillatory  ablative  RM  instability,  presents  a  new  challenge  to  the  imaging 
diagnostics.  With  areal  mass  modulation  amplitude  oscillating  in  time,  the  diagnostics 
that  record  a  limited  number  of  images  incrementally  over  a  relatively  long  period  of  time 
are  difficult  to  apply:  it  is  easy  to  lose  important  information  or  even  completely  miss  the 
effect.  To  extend  the  capabilities  of  our  diagnostic  technique,  we  modified  the  Nike 
imaging  diagnostic  setup  by  adding  a  streak  camera  to  the  system,  which  made  it  possible 
to  analyze  continuous  time  behavior  of  the  x-ray  images. 


Figure  2,  (a)  Orthogonal  x-ray  imaging  diagnostics  implemented  on  Nike  laser.  Face-on  and  side- 
on  images  are  formed  by  two  separate  backlighters  and  spherically  bent  crystals.  For  a  rippled 
double-foil  target  shown  in  the  figure,  the  side-on  streak  records  (left)  shows  the  time  history  of 
collision,  and  the  face-on  streak  record  (right)  -  the  perturbation  evolution,  (b),  (c)  Face-on  streak 
records  of  a  single  30  pm  thick  smooth  target  driven  by  a  4  ns  Nike  pulse  with  a  foot:  (b)  pure 
plastic  target,  (c)  plastic  target  covered  with  a  450  A  gold  layer. 

The  diagnostic  setup  is  presented  in  Fig.  2(a).  Approximately  500  J  is  delivered 
to  one  or  two  silicon  backlighter  targets,  producing  x-rays  that  backlight  the  main  target 
for  about  5  ns.  The  spherically  curved  quartz  crystal  with  the  cut  1011  and  radius  of 
curvature  of  200  mm  selects  the  resonance  line  of  the  He-like  Si  (1.86  keV)  and  projects 
a  monochromatic  image  of  the  target  on  the  slit  of  the  x-ray  streak  camera.  Face-on  and 
side-on  images  are  formed  by  two  separate  backlighters  and  spherically  bent  crystals.  The 
face-on  streak  record  (Fig.  2(a),  right)  shows  the  2D  evolution  of  areal  mass  perturbations 
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in  the  target.  The  monochromatic  x-ray  imaging  makes  it  possible  to  translate  the 
observed  modulation  of  the  optical  thickness  directly  into  the  modulation  of  the  areal 
mass.  Two  examples  of  face-on  streak  records  are  shown  on  Fig.  2(b),  (c).  The  side-on 
streak  records  (Fig  2(a),  left)  show  the  ID  time  history  of  the  CH  foil  acceleration, 
plastic-on-foam  collision,  the  propagation  of  shock  waves  and  the  dynamics  of  the 
ablation  front.  The  streak  records  were  taken  with  a  time  resolution  of  170  ps,  which  is 
sufficient  for  the  -0.5  ns  characteristic  times  of  interest.  A  somewhat  similar  setup  with  a 
thin  slit  as  an  imaging  device  has  been  used  at  NRL  for  the  pioneering  measurements  of 
the  RT  growth  rates  (Whitlock  et  al.  1987). 

To  resolve  the  RM-type  effects  triggered  by  the  well-defined  initial  single-mode 
perturbations  in  the  target  (the  ripples),  it  is  necessary  to  ensure  that  the  driving  radiation 
field  is  sufficiently  smooth,  and  the  non-uniformity  of  the  drive  is  not  an  issue  for  the 
relevant  range  of  perturbation  wavelengths.  This  is  the  case  for  indirectly-driven 
experiments  (Dimonte  &  Remington  1993,  Dimonte  et  al.  1996,  Smitherman  et  al.  1999). 
The  price  paid  for  high  uniformity  of  the  driving  soft  x-rays  in  the  hohlraum  is  its  low 
energy  efficiency:  only  -10%  of  the  total  x-ray  energy  is  delivered  to  the  ablator.  On  the 
other  hand,  if  the  target  is  directly  driven  by  the  laser  radiation,  it  needs  to  be  smoothed, 
which  is  a  general  requirement  for  the  direct-drive  ICF  (Bodner  et  al.  1998,  Lehmberg  & 
Rothenberg  2000).  The  ICF  glass  lasers  use  smoothing  by  spectral  dispersion  (2-D  SSD) 
with  distributed  phase  plate  and  polarization  smoothing  utilizing  birefringent  wedges,  see 
Regan  et  al.  2005  and  references  therein.  For  N  =  10  OMEGA  beams  overlapped  on  the 
target  (Gotchev  et  al.  2006)  the  1.77%  non-uniformity  of  a  single  1  ns  laser  pulse  at  1 

THz  SSD  bandwidth  decreases  by  a  factor  of  ViV  to  anm  =  0.6% .  The  KrF  Nike  laser 
beams  are  smoothed  using  the  induced  spatial  incoherence  (ISI)  technique,  which 
provides  the  best  beam  uniformity  that  has  been  achieved  so  far:  arms  =  1  %  for  a  4-ns 

Nike  beam  at  1  THz  bandwidth  (Deniz  et  al.  1998).  The  RM-type  hydrodynamic 
experiments  on  Nike  (Aglitskiy  et  al.  2001a,  2001b,  2002,  2006,  Weaver  et  al.  2004)  are 
typically  made  with  37-40  beams  overlapping  on  the  target,  which  reduces  arms  below 
0.25%. 

Figures  2(b),  (c)  demonstrate  that  even  such  ultra-uniform  illumination  does  not 
completely  eliminate  the  laser  imprint.  They  present  face-on  streak  records  obtained  for 
30  pm  thick  flat  CH  targets  driven  by  a  4-ns  50  TW/cm  Nike  pulse  with  a  1.6  ns,  8 
TW/cm2  foot.  One  of  these  targets  is  pure  plastic  [Fig.  2(b)],  and  the  front  side  of  the 

O 

other  is  covered  with  a  450  A  layer  of  gold  [Fig.  2(c)],  Both  records  show  no  visible  non¬ 
uniformity  up  to  t  =  1 .5  ns.  Within  this  time  interval,  the  foot  of  the  laser  pulse  pre¬ 
compresses  the  target,  and  its  acceleration  by  the  main  Nike  pulse  starts  at  0.7  ns.  The  RT 
instability  amplifies  the  random  non-uniformity  of  the  laser  irradiation,  and  by  t  =  2  ns 
the  plastic  target  is  already  severely  distorted  by  the  perturbation  growth,  as  evidenced  by 
the  streak  record.  Figure  2(c)  shows  how  the  thin  gold  layer  makes  the  target  more 
imprint-resistant,  as  discussed  by  Obenschain  et  al.  2002:  no  RT  growth  is  observed  till 
the  end  of  the  main  pulse.  To  minimize  the  effect  of  the  random  perturbations  due  to  laser 
imprint,  as  shown  in  this  example,  special  measures  should  be  taken  in  the  RM-type 
hydrodynamic  experiments.  One  needs  either  to  avoid  or  shorten  the  acceleration  phase, 
when  random  laser  perturbations  are  exponentially  amplified  (this  implies  the  use  of 
thicker  targets  or  laser  pulses  without  a  foot)  or  to  use  some  technique  for  imprint 
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mitigation.  In  particular,  in  Section  3  we  discuss  the  effect  of  high-Z  layer  on  the  ablative 
RM  instability. 

To  correctly  model  the  RT  growth,  we  need  to  answer  the  following  important 
question:  How  large  is  the  ablation  front  ripple  amplitude  at  the  moment  when  the 
acceleration  starts?  More  generally,  how  does  this  amplitude  evolve  with  time  during  the 
shock-rarefaction  transit,  before  the  ablation  front  starts  to  accelerate?  Physical 
mechanisms  shaping  this  evolution  are  responsible  for  seeding  the  exponentially  growing 
RT  eigenmodes  at  the  ablation  front.  Their  growth,  in  turn,  limits  the  fusion  energy  gain 
in  the  target. 


-40  -20  0  20  40  -40  -20  0  20  40 

x  (nm)  jt  (urn) 

Figure  3.  RT  seeding  in  a  DT-ice  target  via  outer  (left)  and  inner  (right)  surface  roughness. 
Density  maps  are  shown  after  the  ignitor  shock  launch  (top)  and  as  bum  begins  (bottom). 

This  is  illustrated  by  Fig.  3  showing  density  profiles  for  an  imploding  spherical 
DT  target  obtained  in  a  high-resolution  2-D  simulation  (spherical  harmonics  with  l  from 
1  to  256  are  resolved).  The  simulation  has  been  done  for  the  shock  ignition  direct-drive 
implosion  (Betti  et  al.,  2007),  which  involves  a  fast  re-compression  of  an  expanding 
central  hot  spot  with  a  converging  shock  wave  driven  by  a  short,  powerful  “ignitor”  laser 
pulse  irradiating  the  target  at  peak  compression.  The  two  left  frames  are  simulated  for  the 
perturbations  seeded  by  the  0.49  pm  rms  random  non-uniformity  of  the  outer  target 
surface.  The  right  frames  correspond  to  1  pm  rms  random  non-uniformity  of  the  inner 
target  surface.  Different  perturbation  seeding  and  development  scenarios  affect  final 
compression  and  gain,  which  is  found  to  be  60  in  the  former  case  and  76  in  the  latter.  Of 
course,  separating  these  contributions  to  the  RT  seeding  is  only  possible  in  the 
simulations;  in  the  implosions  of  the  ICF  targets,  all  sources  of  non-uniformity  contribute 
simultaneously.  To  make  sure  that  simulations  describe  these  seeding  mechanisms 
correctly,  we  do  specially  designed  experiments  where  they  can  be  actually  separated. 
These  experiments  provide  data  and  insight  necessary  to  benchmark  the  codes  and  check 
sanity  of  their  predictions. 
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3.  Ablative  RM  instability 

Let  us  turn  on  a  powerful  laser  instantly  upon  a  rippled  surface  of  a  target.  High 
ablative  pressure  launches  a  shock  wave  into  the  target.  If  instead  of  laser  radiation  we 
had  instantly  applied  to  the  target  surface  and  maintained  a  constant  high  pressure,  the 
perturbation  evolution  would  have  been  completely  equivalent  to  the  classical  light-to- 
heavy  RM  instability,  representing  its  limiting  case  A  -*  +1  (Velikovich  et  al.  1998). 

However,  stability  properties  of  an  ablation  front  are  very  different  from  those  of 
a  material  interface.  The  ablation  front  is  stabilized  by  two  distinct  physical  mechanisms. 
One  of  them,  first  discussed  by  Bodner  1974,  is  due  to  the  mass  flow  through  the  ablation 
front.  This  flow  transfers  the  perturbations  from  the  dense  plasma  near  the  ablation  front 
to  the  low-density  plasma  corona.  There  is  also  another,  stronger  stabilizing  effect.  It  was 
discovered  theoretically  by  Sanz  1994,  identified  and  labeled  “the  rocket  effect”  by  Piriz 
et  al.  1997  (see  also  Piriz  2001).  It  emerges  because,  as  first  noted  by  Bodner  1974,  the 
ablation  front  is  an  isotherm.  When  it  is  perturbed,  and  part  of  it  gets  closer  to  the  hot 
corona,  the  temperature  at  this  part  of  the  ablation  front  does  not  increase,  but  the 
temperature  gradient  in  its  vicinity,  VT ,  does.  This,  in  turn,  increases  the  local  heat  flux 
from  the  hot  corona  to  the  ablation  front,  -  /tVT ,  and  hence  the  rate  of  mass  ablation 
from  it,  thereby  increasing  the  ablative  pressure  and  producing  a  restoring  force  that 
pushes  this  part  of  the  ablation  front  back  from  the  laser  absorption  zone.  The  physics  of 
this  “rocket  effect”  is  discussed  in  detail  by  Piriz  et  al.  1997,  Goncharov  1999,  Piriz 
2001. 

With  both  stabilizing  effects  taken  into  account,  the  growth  rate  of  the  ablative  RT 
instability  at  low  acceleration  g  (which  means  high  Froude  number:  Fr  =  v2a  /  gLa  » 1 , 
where  va  =  ml pa  is  the  ablation  velocity,  m  is  the  rate  of  mass  ablation,  pa  is  the 
density  of  the  shock-compressed  target  at  the  ablation  front,  La  is  the  density  gradient 
length  scale  near  the  ablation  front)  is  given  by: 


I  = 


1 


\  1/2 


gk-—k2v2a 


-2  kv„ 


(3) 


(Sanz  1994,  Piriz  et  al.  1997,  Goncharov  et  al.  1996).  Here,  rD  «1  is  the  effective 
blowoff-plasma-to-ablation  front  density  ratio,  a  parameter  of  the  sharp  boundary  model. 
This  parameter  is  a  decreasing  function  of  the  ripple  wavelength  A  =  2jt  I  k  because 
perturbation  eigenmodes  of  longer  wavelengths  extend  farther  into  the  low-density 
corona  plasma;  its  value  is  estimated  using  the  formulas  of  Betti  et  al.  1998a.  Substituting 
into  (3)  g  =  0  (no  acceleration  of  the  ablation  front  during  the  shock/rarefaction  transit), 
we  obtain: 


T  =  ±iQ  -  2kva ,  (4) 

where  Q  =  kva/  . 

The  complex  frequency  (4)  corresponds  to  damped  oscillations  of  the  ablation 
front  in  the  absence  of  acceleration.  The  imaginary  term  in  (4)  describes  stabilization  of 
the  ablation  front  by  the  restoring  force  due  to  the  “rocket  effect”  (Piriz  et  al.  1997, 
Goncharov  1999,  Piriz  2001).  The  real  term  (exponential  decay)  expresses  the 
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stabilization  due  to  the  mass  flow  through  the  ablation  front  (Bodner  1974).  Decaying 
oscillations  of  the  ablation  front  have  been  first  detected  in  simulations  (Velikovich  et  al. 
1998).  The  theory  of  this  effect  was  given,  and  the  term  “ablative  RM  instability”  coined 
by  Goncharov  1999.  Since  the  “ablative  RM  instability”  (that  is,  classical  RM  instability 
stabilized  by  the  ablative  mechanisms)  manifests  itself  in  decaying  oscillations,  they  will 
also  be  referred  below  as  ablative  RM  oscillations. 

The  oscillations  can  be  directly  seen  by  observing  the  time  evolution  of  the  areal 
mass  modulation  amplitude,  which  is  defined  as 

00  +00 

dm  =  d^  p(r  ,6  ,cp)r2  dr  or  dm  =  6  J/?(x,  y,  z)dx  (5) 

0  -00 

for  spherical  and  planar  geometry,  respectively.  Here,  dF  is  the  modulation  amplitude  of 
the  variable  F  with  respect  to  variation  of  the  transverse  coordinates,  (<9,  cp)  or  ( y,  z)  . 

Compared  to  the  displacement  amplitude  of  the  ablation  front,  dxa ,  the  areal  mass 

modulation  amplitude  dm  is  a  much  more  convenient  perturbation  variable  to  be 
estimated  from  the  simulation  results  and  experimental  data.  The  computed  value  of 
dm  is  a  result  of  integration,  which  smoothes  out  numerical  errors,  whereas  determination 
of  dxa  requires  either  numerical  differentiation,  which  is  notoriously  noisy,  or  some  kind 

of  interface  tracking.  Experimentally,  dm  is  evaluated  from  the  modulation  of  the  optical 
thickness,  which  is  directly  observable  with  face-on  x-ray  radiography.  The 
measurements  of  dxa  can  be  made  only  from  side-on  observations,  which  require  a  very 

precise  alignment  and  are  difficult  for  short  perturbation  wavelengths,  -50  pm  or  less. 
However,  the  measured  or  simulated  amplitude  dm  reflects  not  only  the  interfacial 
growth.  Multiple  processes  and  mechanisms  of  perturbation  evolution,  of  which  the 
interfacial  RM-type  growth  might  or  might  not  be  the  dominant  one,  contribute  to  dm .  To 
analyze  the  experimental  data  and  simulation  results  on  dm ,  we  need  to  understand  the 
physics  behind  all  these  contributions. 


Figure  4.  Time  histories  of  the  normalized  areal  mass  modulation  amplitude  (a)  and  its  time 
derivative  (b)  for  the  ablative  RM  instability.  Dotted  lines  in  (b)  show  the  corresponding  results 
for  classical  RM  instability  growth  in  the  limit  A  =  1 .  Shock  velocity  D  =  4x  106  cm/s,  ablation 
velocity  Va  =  105  cm/s,  ripple  wavelength  A  =  30  pm. 

Figure  4  shows  the  oscillations  of  dm  predicted  by  the  small-amplitude  sharp 
boundary  model.  The  total  amplitude  dm  is  shown  in  Fig.  4(a)  as  a  thick  black  line.  The 
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gray  line  shows  the  contribution  to  it  from  the  ablation  front  ripple,  padxa  (where  pa  is 

the  plasma  density  at  the  ablation  front).  The  contribution  to  dm  from  sonic  and  entropy 
perturbations  are  shown  as  thin  solid  line  labeled  S&E,  and  the  contribution  from  the 
rippled  shock  wave  -  as  a  dotted  line.  Figure  4(b)  shows  time  derivatives  of  dm  and 
padxa  as  thick  black  and  gray  lines,  respectively.  To  emphasize  the  difference  between 
the  ablative  and  classical  RM  instability,  thin  lines  in  Fig.  4(b)  show  the  same  variables 
for  the  classical  RM  case  A  =  1  when  the  ablation  effects  are  neglected,  all  other 
parameters  being  the  same. 

This  example  has  been  constructed  for  the  parameters  roughly  corresponding  to 
the  experiments  with  plastic  targets  on  the  Nike  laser:  shock  speed  D  =  4xl06  cm/s, 
shock  compression  ratio  3,  shock  pressure  11  Mbar,  ablation  velocity  va  =105  cm/s, 

Q 

ripple  wavelength  A  =  30  pm,  initial  target  density  /?0=1.07  g/cm  ,  rD=  0.09.  This 

corresponds  to  the  estimated  period  of  the  ablation  front  oscillations  T  =  /  va  =9  ns, 

consistent  with  Fig.  4(a).  The  ablative  RM  rate  (d  /  dt) dm  is  seen  to  oscillate  around  zero, 
whereas  the  classical  RM  growth  rate  tends  to  a  positive  constant  found  from  Wouchuk 
2001.  The  frequency  of  decaying  oscillations  of  the  rippled  shock  wave  transmitted  into 

O  _  1 

the  target  ~  kcs  «  kD  is  ~10  s  ,  much  higher  than  that  of  the  ablation  front  oscillations. 

Amplitude  of  these  oscillations  is  seen  in  Fig.  4(b)  to  be  about  the  same  for  the  cases  of 
classical  and  ablative  RM  instability,  which  is  not  surprising:  oscillations  of  a  rippled 
shock  wave  are  weakly  coupled  to  the  piston  driving  the  shock.  Such  oscillations  have 
been  first  observed  in  the  ICF-relevant  conditions  by  Endo  et  al.  1995  on  GEKKO-XII 
laser.  Both  face-on  x-ray  radiography  and  streaked  optical  self-emission  shock  breakout 
diagnostic  indicated  a  phase  reversal  of  the  rippled  shock  front  when  the  shock 
propagated  ~  0.9 A  into  the  plastic  target.  This  agrees  well  with  our  Fig.  4(a),  where  the 
shock  travel  distance  corresponding  to  the  first  phase  reversal  of  the  rippled  shock  front 
equals  0.85A . 

Figure  4  demonstrates  that  after  the  first  1  ns  the  contributions  both  to  dm  and  to 
(d  /  dt) dm  from  the  ablation  front  dominate.  Observe,  however,  that  the  growth  of  dm  by 
a  factor  of  ~7  shown  in  Fig.  4(a)  is  noticeably  stronger  than  the  corresponding  growth  of 
the  ablation  front  ripple  amplitude,  dxa ,  which  increases  only  by  a  factor  of  2.3.  To 
explain  this  mismatch,  note  that  at  very  early  time 

tin*  P A  a-  (pa  -  Po  )dxs ,  (6) 

where  dxs  is  the  shock  ripple  amplitude,  and  paI  p0=  3  in  this  example  is  the  shock 
compression  ratio.  When  the  shock  wave  is  launched,  dxa  =  dxs  =  dx0 ,  so  that 
dm  =  pAo  =  dm0 .  The  oscillations  of  the  rippled  shock  front,  however,  have  a  much 
higher  frequency  and  decay  much  faster  than  those  of  the  ablation  front.  The  shock  front 
flattens  within  0.6  ns.  The  ablation  front  reaches  its  peak  amplitude  of  dxa  «  2.3xdxlt  at 
about  t  =  1.2  ns.  The  resulting  peak  value  of  dm  is  estimated  from  (6)  by  neglecting  the 
rippled  shock  contribution:  dm  =  padxa  ~  3 Pi)  x  2.3<5r0  =  6.9 dm0 . 
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To  observe  the  ablative  RM  oscillations,  the  laser  pulse  duration  tL  must  be  long 

enough,  the  perturbation  wavelength  A  short  enough,  and  the  target  thickness  L  -  large 
enough,  so  that  the  ablative  pressure  on  the  target  can  be  maintained  for  sufficiently  long 
time  before  it  starts  to  accelerate.  With  the  aid  of  (4),  (6)  one  can  estimate  the 
requirements  on  L  and  tL  that  must  be  met  to  observe  the  first  phase  reversal,  when  both 
the  areal  mass  modulation  amplitude,  dm  ,  and  the  ablation  ripple  amplitude,  dxc , 
change  sign  (Aglitskiy  et  al.  2002): 


L  > 


2  Dr^A 


^  L  ^ 


3  L 
2D 


(7) 


For  the  conditions  of  Fig.  4,  we  obtain  from  (7):  L  >  80  pm,  t,  >  3  ns,  in  agreement  with 
the  figure.  If  the  perturbation  wavelength  is  longer  and  the  laser  pulse  is  shorter,  as  in  the 
experiments  of  Endo  et  al.  1995  (A  =  60  to  100  pm,  L  =  40  to  80  pm,  rL  =2.2  ns, 
I  =  40  TW/cm2,  Al  =  0.53  pm),  the  observation  time  might  not  be  sufficient  to  see  the 
oscillation  of  the  ablation  front.  Indeed,  the  side-on  radiography  used  by  Endo  et  al.  1995 
did  not  detect  any  growth  of  the  ablation  front  ripples.  They  correctly  attributed  the 
changes  in  dm  observed  with  face-on  radiography  to  the  oscillations  of  the  shock  front, 
which  occur  on  a  faster  time  scale,  cf.  Fig.  4. 
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Figure  5.  (a)  Amplitudes  of  the  dominant  Fourier  mode  vs.  time  for  a  93  jam  (black)  and  65  pm 
thick  (gray)  parts  of  the  target;  horizontal  dotted  line  shows  the  initial  peak-to-valley  amplitude, 
(b)  Amplitude  and  phase  of  the  dominant  Fourier  mode  for  a  99  pm  thick  target.  Thick  gray  line 
shows  the  approximate  experimental  uncertainty.  Phase  reversal  of  ^detected  at  about  3.3  ns. 
The  inset  shows  schematic  of  the  experiment. 

The  first  direct  experimental  observations  of  the  ablative  RM  oscillations  have 
been  made  on  the  Nike  laser  at  intensities  about  60  TW/cm  and  reported  by  Aglitskiy  et 
al.  2001a.  They  are  illustrated  by  Fig.  5.  Here  the  Fourier-mode  amplitude  An  is 
normalized  with  respect  to  the  initial  target  density  p{]  and  expressed  in  pm.  The  inset 
illustrates  schematic  of  the  experiment.  Figure  5(a)  presents  the  results  obtained  for  a 
stepped  target  consisting  of  two  halves,  65  pm  and  93  pm  thick.  The  planar  front  surface 
of  this  solid  CH  target  has  been  rippled  at  A  =  45  pm,  peak-to-valley  ripple  initial 
amplitude  3  pm.  Since  the  two  halves  of  this  target  are  driven  by  the  same  laser  pulse,  the 
stepped  target  design  ensures  that  the  ablative  RM  oscillations  should  proceed  identically 
in  both  halves  during  the  shock-rarefaction  transit  time.  Within  the  experimental  error, 
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which  is  larger  for  the  thicker  half  due  to  higher  attenuation  of  the  backlighting  x-rays, 
this  is  indeed  consistent  with  the  observations.  The  peak  value  of  dm  -  5 dm0  in  both 
cases  corresponds  to  the  increase  of  the  ablation  front  ripple  amplitude  from  the  initial 
value  of  dx()  by  the  factor  of  1.5  or  so.  The  main  contribution  to  the  early-time  peak  five¬ 
fold  increase  of  dm  is  due  to  the  flattening  of  the  rippled  shock  front  transmitted  into  the 
target,  as  explained  above. 

The  shock-rarefaction  transit  time  is  shorter  for  the  thinner  half  of  the  stepped 
target.  At  about  2.5  ns  its  acceleration  starts,  switching  the  perturbation  evolution  into  the 
exponential  ablative  RT  growth.  In  the  thicker  half  of  the  target,  the  ablative  RM 
oscillation  continues.  It  approaches  a  phase  reversal  of  dm ,  but  according  to  (7),  the  4  ns 
Nike  laser  pulse  duration  is  not  sufficient  to  observe  it  at  A  =  45  pm.  Such  phase  reversal 
has  been  directly  observed  in  a  thicker  target,  99  pm,  rippled  from  the  front  surface  at 
shorter  wavelength,  A  =  30  pm,  see  Fig.  5(b).  Time  histories  of  both  amplitude  and  phase 
of  the  areal  mass  modulation  both  indicate  phase  reversal  observed  at  about  3  ns,  as  in 
Fig.  6. 

These  direct  observations  of  the  ablative  RM  oscillations  (Aglitskiy  2001a)  and  of 
the  phase  reversal  in  the  target  rippled  from  the  front  side  prior  to  the  start  of  its 
acceleration  (Aglitskiy  2002)  added  confidence  in  our  understanding  of  basic  physics  of 
the  ablative  RM  oscillations.  To  perform  a  more  stringent  test  of  the  theory,  we  need  to 
vary  its  most  significant  parameter,  the  ablation  velocity  va  :  The  theory  predicts  that  the 
efficiency  of  the  “rocket  effect”  stabilizing  mechanism  responsible  for  such  oscillations 
becomes  more  efficient  with  the  increase  of  va .  According  to  the  planar  steady-state 
model  of  the  ablation  front,  the  ablation  velocity  scales  as 


v. a  p:  k 


-1  n-14/15 


(8) 


where  AL  is  the  laser  wavelength  (Manheimer  et  al.  1982).  For  a  given  laser  driver,  the 
straightforward  way  of  increasing  the  ablation  velocity  is  decreasing  the  density  of  the 
target.  Technology  of  the  target  fabrication  makes  it  possible  now  to  do  experiments  with 
low-density  foam  targets  rippled  on  the  front  (Watari  et  al.  2008).  It  is  more  instructive, 
however,  to  see  the  effect  of  ablation  velocity  increase  due  to  density  reduction  on  the 
same  target  as  in  (Aglitskiy  et  al.  2001a,  2002),  in  the  same  shot.  It  can  be  achieved  by 
using  a  target,  a  part  of  which  is  coated  with  a  very  thin  (-1000  A)  layer  of  a  high-Z 
material,  like  gold  or  palladium.  Laser  radiation  absorbed  by  the  high-Z  material  is 
immediately  converted  into  x-rays  and  re-radiated,  heating,  ablating,  and  thereby 
reducing  the  density  of  the  target’s  outer  layers,  until  the  high-Z  layer  expands  and 
becomes  transparent.  The  corresponding  increase  in  va  has  been  demonstrated  to 
strongly  mitigate  the  ablative  RT  instability  (Obenschain  et  al.  2002).  We  have  directly 
observed  its  effect  on  the  ablative  RM  oscillations. 

Figure  6  compares  observed  and  simulated  amplitude  of  the  dominant  Fourier 
mode  vs.  time  for  the  uncoated  and  coated  with  a  800 A  Pd  layer  halves  of  a  40  pm  thick 
plastic  target  rippled  on  the  surface  with  A  =  30  pm  wavelength  and  2.5  pm  amplitude. 
For  the  radiation  transport  in  the  hydro  simulation,  non-LTE  opacity  model  (Klapisch  et 
al.  1998)  was  used.  The  x-ray  heating  resulted  in  a  ~3-fold  decrease  of  density  at  the 
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ablation  front,  which,  according  to  (12)  resulted  in  increasing  the  ablative  velocity  to 
3xl05cm/s. 

The  initial  amplitude  of  the  ripples  is  the  same,  but  their  growth  is  seen  to  be 
more  pronounced  in  the  right,  pure  CH  half  of  the  target.  The  experiment  demonstrates, 
and  simulations  reproduce  lower  amplitude  and  higher  frequency  of  the  oscillations  in  the 
coated  part  of  the  target.  According  to  the  theory,  the  ablation  front  ripple  amplitude  in 
this  part  does  not  grow  at  all;  the  small  observed  growth  of  dm  is  fully  due  to  the 
flattening  of  the  rippled  shock  front.  With  a  higher  oscillation  frequency,  the  theory 
predicts  the  phase  reversal  of  dm  to  occur  at  about  1.5  ns,  much  earlier  than  that  observed 
in  Fig.  6(b).  Due  to  lower  oscillation  amplitude,  the  experiment  did  not  resolve  this  phase 
reversal,  but  the  simulations  confirm  that  it  happens  between  1.5  and  2.5  ns. 


Figure  6.  Observed  (thick  lines)  and  simulated  (thin  lines)  amplitude  of  the  dominant  Fourier 
mode  vs.  time  for  the  uncoated  and  coated  parts  of  the  plastic  target. 

Similar  experiments  with  40  pm  thick  plastic  targets  rippled  at  A  =  20  and  30  pm 
have  been  later  made  on  the  OMEGA  laser  (Gotchev  et  al.  2006),  with  Kirkpatrick-Baez 
microscope  used  as  an  imager  for  face-on  x-ray  radiography.  These  experiments  have 
been  made  at  irradiation  intensity  of  420  TW/cm  ,  much  higher  than  in  (Aglitskiy  2001a, 
2002),  which  made  it  possible  to  detect  the  phase  reversal  of  dm  within  1  ns.  Results  of 
these  accurate  measurements  turned  out  to  be  incompatible  with  modeling  of  laser-target 
interaction  incorporating  flux-limited  Spitzer  thermal  conductivity:  No  choice  of  the  flux 
limiter  value  /  could  make  the  simulations  match  both  the  ID  measurements  of  shock 
breakout  time  and  the  2D  observations  of  the  dm  evolution.  The  matching  had  only  been 
achieved  with  the  aid  of  nonlocal  model  of  electron  thermal  transport  (Goncharov  et  al. 
2006).  The  data  obtained  from  observations  of  the  ablative  RM  oscillations  therefore 
provided  the  first  indication  of  the  nonlocal  character  of  heat  transport  in  laser-driven 
targets. 


4.  Feedout 

Consider  a  planar  shock  wave  breaking  out  at  the  rippled  rear  surface  of  a  target, 
behind  which  there  is  no  low-density  material,  only  vacuum,  or  a  low-density  gas  whose 
presence  does  not  affect  the  hydrodynamics  of  the  shocked  heavy  fluid.  This  corresponds 
to  the  limiting  case  A  =  -1  of  heavy-to-light  classical  RM  instability.  Such  shock-free- 
surface  interaction  does  not  produce  any  material  interface  or  transmitted  shock  wave. 
Only  a  rarefaction  wave  reflected  back  into  the  shocked  target  material  is  produced.  This 
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particular  case  of  a  strong  rarefaction  wave  that  decompresses  to  zero  pressure  is 
sometimes  referred  to  as  an  expansion  wave. 

When  a  planar  or  spherical  shock  wave  arrives  from  the  irradiated  front  or  outer 
surface  of  a  target  to  its  rippled  rear  or  inner  surface,  the  shocked  target  is  turned  into  a 
gas/plasma  (no  material  strength  remains),  and  a  rippled  rarefaction  wave  is  reflected  into 
the  shock-compressed  target.  It  brings  perturbations  to  the  ablation  front,  thereby  feeding 
the  seeds  of  the  subsequent  RT  growth  (cf.  Fig.  3,  right).  This  mechanism  of  RT  seeding 
by  a  rippled  rarefaction  or  expansion  wave  is  called  feedout.  It  had  been  studied  in  detail 
in  planar  geometry  (Betti  et  al.  1998b,  Smitherman  et  al.  1999,  Shigemori  et  al.  2000, 
Velikovich  et  al.  2001,  Aglitskiy  et  al.  2001b,  2002).  Here  is  how  the  feedout  proceeds. 

The  shock  wave  breaking  out  at  the  rear  surface  of  the  ICF  target  is  typically 
strong  enough  for  complete  vaporization  of  the  shocked  material  after  its  decompression. 
The  shocked  material  therefore  expands  as  a  gas,  forming  a  centered 
rarefaction/expansion  wave  (see  Zel’dovich  &  Raizer  2002,  pp.  33-43,  766-770).  The 
head  of  the  expansion  wave  propagates  back  into  the  shock-compressed  material  with 
local  speed  of  sound  cs ,  whose  value  is  uniform  in  the  shocked  target  volume  because  the 
non-uniformity  is  localized  at  its  rear  surface.  Its  tail  expands  at  the  constant  velocity 
vexP  =  2ct  Ky  - 1)  in  the  opposite  direction,  into  the  vacuum,  representing  the  expansion 
front,  where  the  density  and  speed  of  sound  vanish. 

When  the  planar  shock  wave  arrives  to  the  rippled  rear  surface  of  the  target,  it 
first  reaches  the  valleys,  and  later  the  peaks  of  the  ripples.  The  time  delay  between  these 
two  instants  is  2<5x0  /  D  ,  where  2 &c0  is  the  initial  peak-to-valley  ripple  amplitude,  and  D 

is  the  shock  velocity.  The  expansion  starts  from  the  valleys,  while  the  shock  wave 
continues  propagating  through  the  dense  target  material  before  reaching  the  peaks. 
Therefore  the  leading  edge  of  the  rippled  reflected  rarefaction  wave  always  has  the  same 
phase  as  the  initial  surface  ripple.  The  leading  edge  first  reaches  the  planar  front  surface 
of  the  target  where  it  was  initially  thinner,  triggering  the  RT  growth  in  phase  with  the 
initial  distribution  of  the  areal  mass  non-uniformity  dm :  the  initially  thin  and  thick  parts 
of  the  target  grow  thinner  and  thicker,  developing  into  the  bubbles  and  spikes, 
respectively,  at  the  ablation  front,  as  predicted  by  Betti  et  al.  1998b.  Such  development 
had  indeed  been  observed  in  face-on  radiographic  measurements  in  experiments  on 
Gekko-XII  laser  with  25  pm  thick  plastic  targets  rippled  at  the  rear  surface  with  A  =  100 
pm  (Shigemori  et  al.  2000). 

On  the  other  hand,  the  velocity  of  free  expansion  of  the  shocked  plasma  into 
vacuum  typically  exceeds  the  velocity  D  of  the  shock  front  propagating  in  the  same 
direction.  Therefore  while  the  shock  front  propagates  to  the  peaks,  the  expansion  front 
originated  from  the  valleys  is  ahead  of  it.  The  fully  formed  expansion  front  thus  reverses 
its  phase  with  respect  to  the  initial  surface  ripples,  peaks  and  valleys  changing  places. 
Indeed,  such  phase  reversal  of  the  expansion  front  had  been  observed  with  side-on  x-ray 
radiography  by  Shigemori  et  al.  2000. 

If  the  target  thickness  is  not  small  compared  to  the  ripple  wavelength,  then  the 
rippled  rarefaction  wave  itself  affects  the  distribution  of  the  areal  mass  in  the  target. 
When  the  shock  front  breaks  out  at  the  valley,  the  decompression  starts  from  there,  while 
the  high  pressure  is  still  maintained  in  the  shock  wave  propagating  to  the  peaks.  The 
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resulting  lateral  pressure  gradient  starts  driving  the  mass  from  the  peaks  to  the  valleys, 
decreasing  the  pressure  at  the  peak  locations  and  increasing  it  near  the  valley.  The  mass 
flow  continues,  overshooting  the  pressure  equilibrium  situation  and  building  up  a 
reversed  pressure  gradient  (Velikovich  et  al.  2001).  These  sonic  oscillations  of  the  areal 
mass  continue  until  the  leading  edge  of  the  expansion  wave  breaks  out  at  the  ablation 
front,  triggering  the  RT  growth.  The  key  difference  between  such  oscillations  in  shock 
and  rarefaction  waves  is  that  the  restoring  force  is  weaker  in  the  rarefaction  wave,  hence 
the  amplitude  of  the  areal  mass  oscillations  is  much  larger.  For  ideal  gas  EOS,  these 
oscillations  are  described  by  explicit  analytical  formulas  (Velikovich  et  al.  2001,  2005). 
At  late  time  and  for  gas  y  not  too  close  to  unity  but  less  than  2,  these  oscillations  are  well 
described  by  the  asymptotic  formula 


dm 

dm0 


472^-1  +  72^-1)] 

{y-lfQt 


sin  Qt , 


(9) 


where 


Time  (ns)  Distance  (pm)  Time  (ns) 


Figure  7.  (a)  Time  histories  of  the  normalized  areal  mass  modulation  amplitude  in  a  rippled 
expansion  and  shock  waves  for  shock  velocity  D  =  4xl06  cm/s,  ripple  wavelength  A  =  30  jam, 
all  other  parameters  being  the  same  as  in  Fig.  3.3.  The  origin  t  =  0  corresponds  to  the  start  of  the 
laser  pulse  for  the  rippled  shock  wave  and  to  the  moment  of  shock  breakout  at  the  rear  surface  of 
the  target  for  the  rippled  expansion  wave,  (b)  Streak  record  of  a  target  rippled  on  the  rear  side 
with  two  wavelengths:  30  jam  and  45  jam.  (c)  Amplitude  (black)  and  phase  (gray)  of  the 
dominant  Fourier  mode  vs.  time  for  the  45  jam  wavelength.  The  inset  shows  schematic  of  the 
experiment. 

They  are  illustrated  here  by  Fig.  7(a)  plotted  for  y  =  5/3  and  the  same  shock 

velocity  D  =  4xl06  cm/s  as  in  the  example  of  Fig.  5.  For  comparison,  oscillations  of  the 
shock  front  transmitted  into  the  target  in  the  conditions  of  Fig.  5  are  shown  on  the  same 
figure.  The  rarefaction  and  shock  contributions  to  dm  have  been  normalized  differently 
here,  to  facilitate  the  comparison.  For  the  rippled  expansion  wave  dm0  =  p0Sx0,  where 

&0is  the  initial  amplitude  of  the  rear  surface  ripple;  for  the  shock  wave 
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dm0  =  {pa  -  p{)  )dx() ,  where  dxt)  is  the  initial  amplitude  of  the  front  surface  ripple.  The 

shock  front  ripple  is  seen  to  flatten  out  once  the  shock  propagates  one  ripple  wavelength. 
The  rippled  expansion  wave  of  Fig.  7(a)  reverses  the  phase  of  dm  when  its  leading 
propagates  a  quarter  of  the  ripple  wavelength.  Later,  when  the  leading  edge  propagates 
two-thirds  of  the  ripple  wavelength,  the  areal  mass  modulation  reaches  its  first  negative 
peak,  about  3  times  its  initial  amplitude  with  inverted  phase. 

For  these  large  oscillations  of  dm  to  be  observable,  the  ratio  of  target  thickness  L 
to  the  ripple  wavelength  A  should  be  large  enough:  L/ A  >  0.9  to  observe  the  first  phase 
reversal,  LI  A>  2.7  to  observe  the  first  peak  of  dm  with  inverted  phase  (Aglitskiy  et  al. 
2002).  For  the  targets  used  by  Shigemori  et  al.  2000,  LI  A  =  0.25 ,  which  means  that  their 
targets  were  too  thin  to  see  the  dip  in  dm .  The  phase  reversal  of  ^nin  a  rippled 
expansion  wave  was  first  observed  in  the  experiments  with  85  pm  thick  Al  targets  rippled 
on  the  rear  side  at  A  =  50  pm  and  driven  by  the  Nova  hohlraum  (Smitherman  et  al.  1999). 
When  the  head  of  the  rippled  expansion  wave  breaks  out  at  the  ablation  front,  it  starts  the 
target  acceleration  and  triggers  the  RT  growth  of  in  the  dm  in  the  positive  direction,  as 
explained  above.  If  dm  has  already  changed  phase  in  the  rippled  rarefaction  wave  due  to 
the  lateral  mass  flow  in  the  volume  of  the  target,  but  has  not  yet  reached  its  maximum 
with  inverted  phase  at  the  onset  of  acceleration,  then  it  must  change  the  direction  of 
growth,  and  eventually  change  phase  again  due  to  the  areal  mass  redistribution  at  the 
ablation  front.  Such  behavior  had  been  seen  in  the  simulations  of  Smitherman  et  al.  1999 
but  not  resolved  experimentally. 

The  first  direct  observation  of  the  feedout  in  the  regime  LI  A  >  1  where  the  areal 
mass  redistribution  in  the  rippled  expansion  wave  is  important  have  been  reported  by 
Aglitskiy  et  al.  2001b,  2002.  We  observed  the  whole  evolution  of  perturbations  caused  by 
feedout  in  a  planar  target:  large-amplitude  areal  mass  oscillation  in  a  rippled  expansion 
wave,  two  phase  reversals  and  a  subsequent  RT  growth.  Figure  7(b)  shows  the  streak 
record  obtained  for  a  60  pm  thick  plastic  target  rippled  on  the  rear  side  at  A  =  30  pm  and 
45  pm  for  its  left  and  right  half,  respectively.  The  light  and  dark  stripes  indicating, 
respectively,  low  and  high  areal  mass,  are  seen  at  early  time  and  then  they  disappear. 
Later  the  stripes  reappear,  but  the  light  and  dark  ones  change  places,  indicating  a  phase 
reversal,  similar  to  that  shown  in  figure  5(b).  Shortly  after  this,  the  RT  growth  of  dm  in 
the  positive  direction  begins,  and  the  phase  of  dm  is  reversed  again.  Figure  7(c)  shows  the 
time  histories  of  the  peak-to-valley  amplitude  and  phase  of  the  dominant  Fourier  mode, 
normalized  as  above,  for  the  left,  short-wavelength  half  of  the  target.  The  two  phase 
reversals  of  dm  are  clearly  observed. 

5.  Shock  deceleration  of  an  ablation  front 

All  the  experiments  discussed  in  Section  3  refer  to  the  “light-to-heavy”  case  of  the 
ablative  RM  instability,  when  the  shock  wave  originates  at  the  ablation  front  and 
propagates  from  it  into  the  dense  target  plasma.  According  to  the  theory  (Goncharov 
1999,  Goncharov  et  al.  2006),  stability  of  a  non-accelerated  ablation  front  is  an  intrinsic 
property.  Therefore  the  theory  predicts  that  the  ablation  front  should  remain  stable  even 
when  the  shock  wave  arrives  from  the  dense  target  plasma  (the  “heavy-to-light”  case). 
Similarly  to  the  feedout,  the  process  we  are  interested  in  starts  when  a  shock  wave  breaks 
out  at  the  plasma  boundary.  In  this  case,  however,  the  boundary  represents  the  ablation 
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front,  where  a  constant  pressure  is  maintained  by  the  laser,  rather  than  the  free  rear 
surface  at  zero  pressure. 

Our  experiments  were  made  with  targets  consisting  of  planar  plastic  foils  rippled 
on  the  rear  side  (single-mode  ripple  wavelength  A  =  45  pm,  peak-to-valley  amplitude 
2 dx0  =  2  pm)  and  a  planar  layer  of  di vinyl  benzene  foam  (10-20%  of  solid  CH  density) 
separated  from  the  CH  foil  by  a  vacuum  gap,  see  Fig.  2(a).  The  front  side  of  the  plastic 
foil  is  irradiated  by  a  4  ns  Nike  pulse  with  peak  intensity  about  50  TW/cm2.  Laser 
radiation  accelerates  the  CH  foil  across  a  vacuum  gap  until  it  collides  with  a  lower- 
density  foam  layer.  While  the  target  is  accelerated,  a  fast  Rayleigh-Taylor  (RT)  growth  of 
the  single-mode  perturbation  at  the  ablation  front  is  observed.  After  the  collision,  the 
velocity  of  the  ablation  front  is  seen  to  remain  constant,  and  the  RT  growth  is  quenched. 


100  200  300 

Distance  (jim) 


200  400 

Distance  (|im) 


Figure  8.  (a)  Simulated  density  map  and  (b)  side-on  x-ray  streak  record  (b)  for  a  CH  foil  collision 
with  a  foam  layer,  (c)  Side-on  x-ray  streak  record  for  a  directly  irradiated  10%  solid  density  foam 
layer,  (d)  Same  for  a  CH  foil  colliding  with  such  a  layer. 

Figures  8(a),  (b)  demonstrate  the  simulated  and  observed  ID  time  history  of  a 
collision  for  a  target  with  a  30  pm  thick  CH  foil  separated  by  a  120  pm  gap  from  a  200 
pm  thick,  10%  solid  CH  density  foam  layer.  Figure  8(a)  is  a  simulated  density  map  on 
the  (x,t)  plane.  The  laser  irradiation  drives  a  strong  shock  wave  into  the  CH  foil, 
compresses  it  and  sets  it  into  the  motion  at  constant  mass  velocity  that  is  associated  with 
the  shock  wave.  Before  the  start  of  its  acceleration,  the  trajectory  of  the  ablation  front  on 
the  (x,t)  plane  is  straight.  The  compressed  CH  foil  plasma  starts  to  accelerate  following 
the  break-out  of  this  shock  at  its  rear  surface,  when  the  rarefaction  wave  reflected  from  it 
reaches  the  ablation  front  ( t « 1 .4  ns).  The  acceleration  makes  the  trajectory  of  the 
ablation  front  on  the  (x,t)  plane  convex.  Then  at  t  ®  3  ns  the  foil  accelerated  to  ~107  cm/s 
collides  with  the  foam  layer.  Simulations  indicate  that  the  collision  briefly  produces  a 
pressure  of  ~30  Mbar,  which  is  3.5  times  greater  than  the  pressure  ~8  Mbar  maintained 
by  the  laser  at  the  ablation  front.  This  pressure  pulse  generates  two  strong  shock  waves. 
One  of  them  propagates  forward  into  the  thick  foam  layer,  while  the  other  propagates 
back  into  the  thin,  compressed  and  partly  ablated  CH  foil  plasma.  The  ablation  front  is  re¬ 
shocked  from  the  dense  plasma  side.  The  ablation  front  rapidly  decompresses  from  the 
high  shock  pressure  to  ~8  Mbar.  The  rarefaction  wave  that  emerges  from  the  ablation 
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front  at  the  instant  of  its  re-shock  carries  the  reduced  pressure  back  to  the  shock  wave 
propagating  into  the  foam,  rapidly  slowing  it  down. 

This  shock  is  then  fully  determined  by  the  ablative  pressure  maintained  by  the 
laser  ablation  and  the  foam  density,  as  evidenced  by  the  shape  of  this  shock’s  ( x,t ) 
trajectory  on  Fig.  8(a),  which  is  concave  for  a  short  time  interval,  and  then  a  straight  line 
after  t  =  3.2  ns.  The  re-shock  impulsively  changes  the  velocity  of  the  ablation  front, 
which  then  stays  constant,  and  its  ( x,t )  trajectory  becomes  straight  again  starting  from 
r  =  3.1  ns.  Figure  8(c),  (d)  compare  velocities  of  the  shocks  in  10%  solid  density  foam 
produced  by  direct  irradiation  of  the  foam  with  Nike  beams  (a)  and  by  collision  of  the 
laser-accelerated  CH  foil  with  the  foam  layer  (b).  The  slope  of  the  dashed  lines  exactly 
corresponds  to  the  shock  velocity  in  the  latter  case,  averaged  from  the  instant  of  collision 
to  the  observed  shock  breakout  at  the  rear  surface,  and  it  is  very  close  to  the  observed 
velocity  of  the  shock  wave  driven  into  the  foam  directly  by  the  laser.  The  constant 
velocities  of  the  shock  front  in  the  foam  and  of  the  ablation  fronts  at  this  stage  are  fully 
determined  by  the  ablative  pressure  maintained  by  the  laser  and  the  foam  density. 

Observed  and  simulated  evolution  of  the  dominant  Fourier  mode  of  the  areal  mass 
modulation  is  shown  in  Fig  9(a).  The  thick  lines  with  approximate  error  bars  present  the 
experimental  data  from  the  face-on  images.  The  pink  line  corresponds  to  the  same 
conditions  as  Fig.  8(a),  (b),  the  light  blue  line  to  a  higher  foam  density:  10%  and  20%  of 
solid,  respectively.  The  signal  is  normalized  with  respect  to  its  initial  value  dm0  =  p0dx0 . 
The  small  time  shift  between  these  two  curves  prior  to  the  collision  is  due  to  a 
combination  of  difference  in  the  actual  foil  thickness  and  the  timing  relative  to  the  start  of 
the  laser  pulse.  Thin  red  and  blue  lines  show  the  simulation  results  obtained  in  two- 
dimensional  simulations  for  the  same  conditions  as  pink  and  light  blue  lines,  respectively, 
with  plasma  radiation  taken  into  account. 


Figure  9.  (a)  Observed  (thick,  with  approximate  error  bars)  and  simulated  (thin  lines)  time  history 
of  dm  in  the  target.  Arrow  indicates  the  instant  of  re-shock  on  simulated  lines,  (b)  Simulated  time 
history  of  the  dm  in  a  CH  foil  freely  accelerated  (green  line)  and  colliding  with  a  foam  layer  (red 
and  blue  lines),  (c)  Observed  time  history  of  dm  in  the  CH  foil  rippled  on  the  rear  side  (green  line 
with  approximate  error  bars)  and  in  such  foil  colliding  with  a  foam  layer  (light  blue  line). 

Figure  9(a)  shows  that  the  early-time  minimum  and  two  phase  reversals 
characteristic  of  the  feedout.  They  are  not  well-resolved  here  because  our  CH  foil  is 
thinner  than  those  used  in  the  feedout  experiments,  see  Section  4.  Our  observations  are 
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nevertheless  consistent  with  the  presence  of  this  minimum  and  with  the  corresponding 
delay  in  the  onset  of  the  RT  growth.  It  is  thereby  seeded  by  the  feedout  mechanism,  see 
Section  4. 

The  RT  growth  is  quenched  when  the  collision  occurs  and  the  ablation  front  is  re¬ 
shocked.  A  re-shock  of  a  material  interface  would  (re)start  a  classical  RM  growth  of  its 
perturbations.  Since  the  linear  RM  growth  rate  is  proportional  to  the  initial  ripple 
amplitude,  a  “heavy-to-light”  re-shock  of  a  material  interface  after  a  substantial  RT  or 
RM  amplification  of  its  non-uniformities  would  give  rise  to  a  rapid  growth  of  the  areal 
mass  modulation  amplitude  in  the  negative  direction,  implying  a  quick  reversal  of  its 
phase.  A  similar  increase  of  the  RM  growth  rate  at  a  re-shocked  interface  has  been 
observed  in  shock-tube  experiments  (Leinov  et  al.  2008).  Our  results  presented  in  Fig. 
9(a)  do  not  demonstrate  such  a  behavior.  Instead,  the  growth  of  the  areal  mass 
modulation  amplitude  is  seen  to  end  with  the  re-shock  and  to  be  followed  by  decaying 
oscillations  around  a  nonzero  average  value.  This  effect  is  robust  and  caused  entirely  by 
the  collision. 

As  demonstrated  in  Figs.  9(b),  (c),  the  quenching  is  not  due  to  nonlinear 
saturation  or  any  other  mechanism  that  would  slow  down  the  RT  growth  in  the  absence 
of  the  collision.  Without  the  collision,  the  simulated  RT  growth  (green  line  in  Fig.  9(b)) 
would  continue  until  the  target  bums  through,  or  is  broken  up  by  the  instability,  or  the 
laser  pulse  ends.  The  same  is  observed  in  the  experiment,  Fig.  9(c).  Here  the  collision 
occurs  at  3.5  ns,  later  than  in  Figs.  9(a),  (b),  due  to  lower  laser  intensity  and  wider  gap 
between  the  foil  and  the  foam,  but  still  well  within  the  4  ns  Nike  pulse. 

Decaying  oscillations  of  a  non-accelerated  ablation  front  seen  in  Figs.  9(a),  (b)  are 
not  the  low-frequency  ablative  RM  oscillations  (cf.  Figs.  4,  5).  These  are  oscillations 
around  a  non-zero  average.  The  ablation  front  ripples,  which  provide  the  largest 
contribution  to  the  observed  areal  mass  modulation  amplitude  dm ,  do  not  change  phase. 
On  the  other  hand,  the  simulated  and  observed  oscillation  period,  ~1  ns,  is  close  to  the 
sound  wave  period  X  /  cs ,  where  A  =  45  pm,  and  the  speed  of  sound  cs  in  the  re-shocked 

foil  plasma  is  estimated  from  the  simulations  to  be  ~5xl06cm/s.  The  observed  time 
history  of  An  therefore  involves  fast  areal  mass  oscillations  in  the  rippled  rarefaction 
wave  that  emerges  from  the  ablation  front  after  the  re-shock  (as  in  the  case  of  feedout,  cf. 
Fig.  7(a)),  superimposed  upon  the  slow,  also  oscillatory,  evolution  due  to  the  ablative  RM 
instability.  Indeed,  all  our  simulation  results  show  that  the  changes  in  dm  after  the 
collision  accumulate  in  the  volume  of  the  plasma  reached  by  the  rarefaction  wave. 

Therefore,  our  observations  prove  that  the  ablation  front  becomes  stable  as  soon 
as  its  acceleration  stops,  in  agreement  with  the  theory. 

6.  Impulsive  Loading 

The  process  called  impulsive  loading  can  be  idealized  as  a  fast  release  of  a  finite 
energy  in  a  thin  layer  near  the  surface  x  =  0  of  a  half-space  x  <  0  filled  with  a  uniform 
material  at  rest.  Impulsive  loading  had  been  studied  since  the  1950s  (Hafele  1955, 
Zel’dovich  &  Raizer  2002  pp.  820-849  and  references  therein),  when  it  was  visualized  as 
an  explosion  of  a  thin  surface  layer  of  explosive,  or  a  concentrated  impact  of  a  thin,  light 
plate  carrying  a  finite  kinetic  energy.  A  short  laser  pulse  depositing  a  finite  thermal 
energy  in  a  thin  surface  layer  of  a  target  before  being  turned  off  also  can  act  as  an 
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impulsive  load.  Impulsive  loading  sends  a  strong  shock  wave  into  the  target,  which  is 
immediately  followed  by  an  expansion  wave  that  gradually  decelerates  it.  In  a  sense,  the 
impulsive  loading  is  similar  to  a  planar  blast  wave  (Sedov  1959),  and  it  also  produces  a 
self-similar  flow  of  the  shocked  material.  The  difference  is  that  the  impulsive  loading 
involves  free  expansion  of  this  material  to  zero  pressure  (in  a  blast  wave,  the  minimum 
pressure  in  a  rarefaction  wave  following  the  shock  wave  is  finite),  which  implies  a 
different  kind  of  self-similarity,  see  the  discussion  by  Zel’dovich  &  Raizer  2002,  Chapter 
XII. 

Renewed  interest  to  the  impulsive  loading  in  the  area  of  direct-drive  laser  fusion 
emerged  several  years  ago.  It  is  due  to  the  recent  proposals  of  using  a  short  laser  pulse 
called  a  spike  or  a  picket  prior  to  the  main  drive  pulse  in  order  to  make  the  target  more 
resistant  to  laser  imprint  (Metzler  et  al.  2002,  2003)  and  to  shape  the  target  adiabat 
(Goncharov  et  al.  2003,  Anderson  &  Betti  2003,  Betti  et  al.  2005).  The  stabilizing  effects 
of  the  impulsive  loading  are  caused  by  the  density  profile  in  the  expansion  wave  that 
follows  the  decaying  shock  wave.  Since  the  density  increases  into  the  target,  the  shock 
wave  driven  by  the  main  pulse  slows  down  as  it  propagates,  and  so  does  the  ablation 
front  pushing  it.  Deceleration  stabilizes  the  ablation  front  much  more  strongly  than  the 
“rocket  effect”  responsible  for  the  ablative  RM  oscillations,  which  results  in  effective 
mitigation  of  perturbation  growth  due  to  the  non-uniformities  of  both  laser  irradiation  and 
target  surface  (Metzler  et  al.  1999).  Moreover,  as  the  shock  wave  driven  by  the  main 
pulse  slows  down,  the  density  of  the  shocked  target  material  increases  and  its  temperature 
decreases.  We  end  up  with  a  “tailored  adiabat”  of  the  target,  where  adiabat  is  defined  as 
a  =  pi pF  p~5'3 ,  p  and  pF being,  respectively,  the  pressure  in  the  target  and  the 
pressure  of  degenerate  electron  Fermi-gas  at  the  given  density  p.  Lower  density  in  the 
outer  layers  of  the  target  ensures  a  higher  ablation  velocity  va ,  hence,  lower  rate  of  the 
RT  growth,  cf.  Eq.  (3),  and  thereby  better  stability  of  implosion.  On  the  other  hand,  high 
adiabat  a  of  the  fuel  makes  it  difficult  to  compress.  The  kinetic  energy  of  the  target 
requires  for  its  ignition  scales  as  a  power  of  the  in-flight  adiabat  aif :  Eig  oc  a-f ,  where 

the  high  power  k  is  estimated  to  be  between  1.8  (Herrmann  et  al.  2001,  Kemp  et  al. 
2001)  to  2.4  (Betti  et  al.  2002).  Using  the  impulsive  loading  to  tailor  the  target  adiabat, 
we  can  achieve  stable  implosion  and  high  compression  of  the  fuel  simultaneously. 


Figure  10.  (a)  Time  evolution  of  a  density  profile  in  a  foam  layer  irradiated  by  a  0.325  ns  long 
pulse.  The  profiles  are  shown  with  a  1  ns  time  interval,  (b)  Simulated  amplitudes  of  the  dominant 
Fourier  mode  vs.  normalized  shock  travel  distance  for  an  impulsively  loaded  rippled  target,  (c) 
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Observed  dominant  Fourier  mode  vs.  time  for  a  rippled  plastic  target  loaded  with  a  0.3  ns  long 
pulse.  The  inset  shows  schematic  of  the  experiment. 

Figure  10(a)  illustrates  the  ID  flow  produced  by  the  impulsive  loading  in  a  planar 
geometry.  It  shows  a  sequence  of  simulated  density  profiles  resulting  from  the  irradiation 
of  a  CH  foam  layer  (density  0.39  g/cm3)  with  a  0.325  ns  long  KrF  pulse  at  4.1  TW/cm2. 
A  graded  density  profile  produced  by  the  shock-expansion  flow  expands  with  time  in 
approximately  self-similar  way.  After  being  shocked,  each  particle  first  moves  in  the 
direction  of  shock  propagation.  Later,  the  negative  pressure  gradient  slows  it  down  and 
accelerates  it  in  the  opposite  direction. 

Both  rippled  shock  and  expansion  waves  tend  to  oscillate  as  they  propagate  at 
constant  velocity  into  an  unperturbed  fluid.  Oscillations  of  the  expansion  wave  are 
stronger,  cf.  Fig.  7(a).  Evolution  of  perturbations  in  the  shock-expansion  flow  generated 
by  impulsive  loading  was  first  studied  theoretically  by  Velikovich  et  al.  2003.  It  was 
found  that  the  coupled  shock  and  expansion  waves  produce  much  stronger  oscillations 
than  either  a  shock  wave  or  an  expansion  wave  separately.  This  happens  because  the 
shock  wave  decelerates,  weakening  the  restoring  force  that  acts  upon  the  shock  ripples. 
As  the  decelerating  shock  wave  propagates  through  the  target,  its  modulation  amplitude 
grows,  exceeding  the  initial  amplitude  by  a  factor  of  2  or  more.  The  oscillating  areal  mass 
modulation  amplitude  dm  reaches  the  peak  values  that  exceed  its  initial  value 
dm0  =  p0dx0  by  a  factor  of  5-7  or  more  and  reverses  its  phase  several  times  after  the  laser 
pulse  is  over. 

This  behavior  is  illustrated  by  Fig.  10(b).  Here  the  areal  mass  modulation 
amplitude  is  shown  vs.  the  normalized  distance  traveled  by  the  shock  wave,  x/A,  for 
ideal-gas  equations  of  state  with  y  varied  from  6/5  to  5/3  and  CALEOS  model  equation  of 
state  of  plastic  foam.  The  oscillations  are  seen  to  be  the  strongest  for  the  most 
compressible  EOS,  y  =  6/5,  which  reveals  the  physical  similarity  of  this  strong  effect 
with  Vishniac’s  instability  of  a  blast  wave  (Vishniac  1983,  Ryu  &  Vishniac  1987).  As 
explained  by  Velikovich  et  al.  2005,  the  Vishniac’s  instability  is  also  driven  by  the 
instability  of  a  planar  expansion  wave  at  low  y. 

Figure  10(c)  shows  the  first  experimental  observation  of  such  oscillations  reported 
by  Weaver  et  al.  2004.  A  0.3  ns  long,  7.7  TW/cm  KrF  laser  spike  irradiated  the  target  3 
ns  prior  to  the  arrival  of  the  main  pulse  at  t  =  0  ,  rapidly  depositing  2.3  kJ/cm  of  energy 
into  a  40  pm  thick  plastic  foil  rippled  on  the  front  surface  with  A  =  30  pm,  dx(j  =  2.5  pm. 
About  t  =  0  the  shock  wave  driven  by  the  spike  breaks  out  at  the  rear  surface  of  the 
target,  which  corresponds  to  kxs  =  8.37  in  Fig.  10(a).  Between  t  =  -2.7  ns  and  t  =  0  the 

laser  is  off,  and  the  evolution  of  dm  seen  in  Fig.  10(c)  are  represents  free  oscillations  of 
the  shock-expansion  flow  described  above.  We  see  the  expected  strong  oscillations  of  the 
areal  mass  and  its  phase  reversal  at  about  the  same  time  as  predicted. 

Detailed  experimental  study  of  the  strong  shock  wave  and  areal  mass  oscillations 
associated  with  impulsive  loading  of  planar  targets  is  still  ahead.  It  will  provide  valuable 
data  for  benchmarking  hydro  codes  and  testing  the  EOS  approximations  in  the  Mbar 
pressure  range.  The  observable  oscillations  take  place  when  the  laser  is  off  and  are  fully 
determined  by  the  integral  characteristics  of  the  short  pulse:  the  energy  deposited  near  the 
target  surface  and  the  preheat  of  the  target  material  by  x-rays  or  fast  electrons,  if  any. 
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Therefore,  study  of  these  oscillations  can  be  largely  separated  from  the  complicated 
physics  of  laser-plasma  interaction. 


7.  Conclusions 

A  considerable  progress  had  been  made  in  the  last  decade  in  observation  and 
understanding  the  RM-type  hydrodynamics  processes.  An  extensive  data  base  has  been 
accumulated  for  benchmarking  hydrocodes  and  testing  theoretical  models.  Agreement 
between  the  experimental  data  and  simulation  result  in  most  cases  is  reasonably  good, 
which  indicates  that  the  basic  physics  underlying  the  RM-type  phenomena  is  understood 
correctly.  A  decisive  role  in  this  progress  has  been  played  by  the  monochromatic  x-ray 
imaging,  which  proved  to  be  a  powerful  tool  for  observing  dynamics  of  ICF  targets  and 
probing  the  lateral  mass  redistribution  due  to  hydro  instabilities,  waves  and  other 
transient  processes. 

Of  course,  much  more  needs  to  be  done  to  fine-tune  the  details  and  to  make  sure 
that  all  the  relevant  physics,  such  as  equation  of  state  and  radiation  transfer/preheat,  is 
properly  accounted  for.  Some  of  the  theoretically  predicted  RM-type  processes  are  yet  to 
be  observed,  for  example,  the  instability  on  the  material  interface  that  is  dynamically 
formed  at  collision  of  two  planar  foils,  when  at  least  one  of  the  colliding  foils  is  rippled 
(Velikovich  et  al.  2000).  In  addition  to  the  existing  ICF  laser  facilities,  the  new  10-kJ 
class  lasers,  such  as  LIL,  Shenguang-III  and  Luch,  open  the  opportunities  for  such 
experiments  to  be  carried  out  in  France,  China  and  Russia.  This  research  will  be  an 
important  contribution  to  the  efforts  aimed  at  achieving  ignition,  fusion  energy  gain,  and 
eventually,  fusion  energy  production,  on  NIF  in  the  U.  S.,  LMJ  in  France  and  next- 
generation  laser  facilities  of  the  future. 
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